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ABSTRACT 

Using the fluid-magnetofluid formalism, we obtain axisymmetric stability criteria for a 
composite disc system consisting of stellar and gaseous magnetized singular isothermal 
discs (MSIDs). Both (M)SIDs are presumed to be razor-thin and are gravitationally 
coupled in a self-consistent axisymmetric background equilibrium with power-law sur- 
face mass densities and flat rotation curves. The gaseous MSID is embedded with a 
non- force- free coplanar azimuthal magnetic field Bg{r) of radial scaling r~^/^. Lou 
& Zou recently constructed exact global stationary configurations for both axisym- 
metric and nonaxisymmetric coplanar magnetohydrodynamic (MHD) perturbations 
in such a composite MSID system and proposed the MHD Dg— criteria for axisym- 
metric stability by the hydrodynamic analogy. In a different perspective, we derive and 
analyze here the time-dependent WKBJ dispersion relation in the low-frequency and 
tight-winding regime to examine axisymmetric stability properties. By introducing a 
rotational Mach number Dg for the ratio of the stellar rotation speed Vg to the stellar 
velocity dispersion Os, one readily determines the stable range of numerically to 
establish the Dg— criteria for axisymmetric MSID stability. Those MSID systems ro- 
tating either too fast (ring fragmentation) or too slow (Jeans collapse) are unstable. 
The stable range of depends on three dimensionless parameters: the ratio A for the 
Alfven speed to the sound speed in the gaseous MSID, the ratio /3 for the square of the 
stellar velocity dispersion to the gas sound speed and the ratio 6 for the surface mass 
densities of the two (M)SIDs. Our WKBJ results of (M)SID instability provide phys- 
ically compelling explanations for the stationary analysis of Lou & Zou. We further 
introduce an effective MHD Q parameter for a composite MSID system and compare 
with the earlier work of Elmegreen, Jog and Shen & Lou. As expected, an axisym- 
metric dark matter halo enhances the stability against axisymmetric disturbances in 
a composite partial MSID system. In terms of the global star formation rate in a disc 
galaxy system, it would appear physically more sensible to examine the MHD Qm 
stability criterion against galactic observations. Relevance to large-scale structures in 
disc galaxies are also discussed. 

Key words: MHD waves — ISM: magnetic fields — galaxies: kinematics and dynam- 
ics — galaxies: spiral — star: formation — galaxies: structure. 



1 INTRODUCTION 

In contexts of galactic structures, disc stabilities and global 
star formation rates in spiral galaxies, we derive here in- 
stability criteria for axisymmetric coplanar magnetohydro- 
dynamic (MHD) perturbations in a composite disc system 



with an azimuthal magnetic field in the gas disc component 
and establish a generalized definition of an effective MHD 
Qm parameter appropriate to such a magnetized gravita- 
tional system. Formulated as such, this is an idealized and 
limited theoretical MHD disc problem yet with several key 
conceptual elements included. The simple physical rationale 
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is that zonal regions of higher gas density and magnetic field 
are vulnerable to active formation of massive stars with var- 
ious scales involved differing by many orders of magnitudes. 
On the same ground, non-axisymmetric stability criteria are 
equally important but are more challenging to establish (see 
e.g. Shu ot al. 2000 for relevant issues). Over four decades, 
important development have been made for instability crite- 
ria relevant to galactic disc dynamics (see Lin 1987, Binney 
& Tremaine 1987 and Berlin & Lin 1996 and extensive refer- 
ences therein) . The original studies of axisymmetric instabil- 
ities were conducted by Safronov (1960) and Toomre (1964) 
who introduced the dimensionless Q parameter for the local 
stability (i.e. Q > 1) against axisymmetric ring-like distur- 
bances. For disc galaxies, it would be more realistic and sen- 
sible to investigate a composite disc system consisting of a 
stellar disc, a magnetized gas disc and a massive dark mat- 
ter halo. There have been extensive theoretical studies on 
this type of composite two-component disc systems of var- 
ious sub-combinations (Lin & Shu 1966, 1968; Kato 1972; 
Jog & Solomon 1984a; Bertin & Romeo 1988; Romeo 1992; 
Elmegreen 1995; Jog 1996; Lou & Fan 1998b, 2000a, b; Lou 
& Shen 2003; Shen & Lou 2003, 2004a, b; Lou & Zou 2004; 
Lou & Wu 2005). In particular, there have been several stud- 
ies trying to define a proper effective Q parameter for local 
axisymmetric instability criterion in a composite disc sys- 
tem (Elmegreen 1995; Jog 1996; Lou & Fan 1998b, 2000a, 
b; Shen & Lou 2003, 2004a, b). From different perspectives, 
these analyses offer insights for instability properties of a 
composite disc system and provide a theoretical basis for 
understanding the large-scale dynamics in such a system 
(e.g. Lou & Fan 2000a, b; Lou & Shen 2003; Shen & Lou 
2003, 2004a, b; Lou & Zou 2004; Lou & Wu 2005). 

The main motivation here is to explore basic proper- 
ties the MSID model in a composite system and obtain con- 
ceptual insights for astrophysical applications in magnetized 
spiral galaxies and in estimating global star formation rates 
in disc galaxies (e.g., Lou & Bian 2005). In theoretical stud- 
ies of modeling disc galaxies, the class of SID models has 
a distinguished history since the pioneering work of Mestel 
(1963) (Zang 1976; Toomre 1977; Lemos, Kalnajs & Lynden- 
Bell 1991; Lynden-Bell & Lemos 1993; Syer & Tremaine 
1996; Goodman & Evans 1999; Charkrabarti, Laughlin & 
Shu 2003) and has gained considerable attention and inter- 
ests recently by considering a composite disc system and by 
incorporating effects of magnetic field (Shu et al. 2000; Lou 
2002; Lou & Fan 2002; Lou & Shen 2003; Shen & Lou 2003, 
2004a, b; Shen, Liu & Lou 2004; Lou & Zou 2004; Lou & Wu 
2005). Specifically, Shu et al. (2000; sec also Galli et al. 2001) 
studied global stationary (i.e., zero pattern speed) perturba- 
tion configurations in an isopedically magnetized SID with- 
out invoking the usual WKBJ or tight-winding approxima- 
tion. They obtained exact global solutions for both aligned 
and unaligned axisymmetric and non-axisymmetric logarith- 
mic spiral configurations and interpreted the axisymmetric 
solution for perturbations with radial propagations as de- 
marcating the boundaries between the stable and unstable 
regimes. By these axisymmetric instabilities, a SID with a 



sufficiently slow rotation speed would Jeans collapse induced 
by perturbations of larger radial scales, while a SID with a 
sufficiently fast rotation speed may suffer the ring fragmen- 
tation instability induced by perturbations of smaller radial 
scales (see fig. 2 of Shu et al. 2000). By introducing a ro- 
tational Mach number D, defined as the ratio of the SID 
rotation speed V to the isothermal sound speed a, the crit- 
ical values of the highest and lowest D for an axisymmetric 
stability can bo determined directly from the marginal sta- 
bility curve. To support their physical interpretations, they 
invoked the well-known Toomre Q parameter and found that 
the highest D, namely the minimum of the ring fragmenta- 
tion curve, corresponds to a Q value very close to unity, 
thus heuristically suggesting the correspondence between 
the 75— criterion and the criterion. 

Different from yet complementary to the analysis of Shu 
et al. (2000) on a single isopedically magnetized disc, Lou 
(2002) studied global coplanar MHD perturbations in a sin- 
gle MSID embedded with an azimuthal magnetic field and 
revealed that the minimum of the MHD ring fragmenta- 
tion curve in this MSID model is tightly associated with 
the generalized MHD Qm parameter originally introduced 
by Lou & Fan (1998a) in developing the galactic MHD den- 
sity wave theory (Fan & Lou 1996). For a composite system 
of two coupled unmagnetized hydrodynamic SIDs, Lou & 
Shen (2003) constructed stationary global perturbation con- 
figurations and Shen & Lou (2003) suggested a straightfor- 
ward D— criterion for the axisymmetric ring fragmentation 
instability in such a system on the basis of a low-frequency 
WKBJ analysis; they revealed that the minimum of the ring 
fragmentation in their composite SID model is again closely 
related to a proper effective Q parameter (Elmegreen 1995; 
Jog 1996). Furthermore, for a composite SID system with 
an isopedically magnetized gaseous SID and a stellar SID 
in the fiuid description, Lou & Wu (2005) have constructed 
global stationary MHD perturbation structures and exam- 
ined stability properties to anticipate a similar D— criterion 
in parallel to the case of Shen & Lou (2003). Meanwhile, 
Shen & Lou (2004b) have further generalized both work of 
Syer & Tremaine (1996) and Lou & Shen (2003) to the situ- 
ation of a composite system for two gravitationally coupled 
scale-free discs; they also studied the axisymmetric stability 
of such a composite system in terms of the marginal stability 
curves and proposed a D^— criterion by the analogy of Shen 
& Lou (2003). 

We have recently examined two-dimensional coplanar 
MHD perturbations in a composite system consisting of a 
stellar SID and a gaseous MSID. Both SIDs arc expediently 
approximated as razor-thin and the gas disc is embedded 
with a non-force-free coplanar magnetic field (Lou & Zou 
2004). In this fluid-magnetofluid MSID model approxima- 
tion, we obtained exact global stationary MHD solutions for 
aligned and unaligned logarithmic spiral perturbation con- 
figurations in such a composite MSID system, expressed in 
terms of the stellar rotational Mach number Dg. In refer- 
ence to the results of a single SID (Shu et al. 2000; Galh 
et al. 2001; Lou 2002), it would be natural to suggest that 
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the stationary axisymmetric solutions with radial propaga- 
tions give rise to marginal stability curves (see Figure |5| in 
this paper later and Lou & Zou 2004). In comparison with 
the single SID case, the stable range of Dj. is reduced as a 
result of the mutual gravitational coupling between the two 
SIDs. However, in comparison with the case of a composite 
unmagnetized SID system (Lou & Shen 2003; Shen & Lou 
2003), the stable range of expands considerably due to 
the presence of a coplanar magnetic field. 

To confirm heuristic arguments for the above analogy 
and our intuitive physical interpretations, we conduct in this 
paper a low-frequency time-dependent stability analysis in 
the WKBJ or tight-winding approximation for the compos- 
ite MSID system (Shen & Lou 2003, 2004a; Lou & Zou 2004). 
We shall demonstrate unambiguously the validity of demar- 
cating the stable and unstable regimes by the stellar rota- 
tional Mach number Ds- To place our analysis in proper 
contexts, we also discuss specifically how the two effective 
Q parameters of Elmegreen (1995) and of Jog (1996) are re- 
lated to our Ds parameter when other relevant parameters 
are specified, and show that the two effective Q parame- 
ters are pertinent to the ring fragmentation instability in a 
composite MSID system. 

In Section 2, we derive the time-dependent disper- 
sion relation using the WKBJ approximation for pertur- 
bations in a composite MSID system and introduce a few 
key dimensionless parameters. In Section 3, we present the 
Ds— criterion and two effective Q parameters for a compos- 
ite MSID system being stable against arbitrary axisymmet- 
ric perturbations with radial propagations. Main results and 
discussions are summarized in Section 4. 



2 A FLUID-MAGNETOFLUID FORMALISM 

We consider below a composite system consisting of two 
gravitationally coupled (M)SIDs with one of the SIDs being 
magnetized and thus referred to as MSID. For physical vari- 
ables, we use either superscript or subscript s to indicate 
an association with the stellar SID and either superscript 
or subscript g to indicate an association with the gaseous 
MSID. The stellar SID and the gaseous MSID can have con- 
stant yet different rotational speeds 14 and Vg (related to 
the phenomenon of asymmetric drift in the galactic context) ; 
we thus write the background angular rotation speeds fis of 
the stellar SID and f2g of the gaseous MSID as 

0.S = Vs/r = asDs/r (1) 

and 

fig = VJr = agDJr , (2) 

separately, where Us and ag are the constant velocity dis- 
persion of the stellar SID and the isothermal sound speed 
of the gaseous MSID, respectively; Ds and Dg are the cor- 
responding rotational Mach numbers. The relevant epicyclic 
frequencies in terms of f2s and Q.g are given by 



and 



^-^(r^n.) = 2J7? 
r dr 



(3) 



r dr 



(4) 



respectively. Similar to a single MSID, we take the back- 
ground azimuthal magnetic field to be in the form of 

Be{r)^Tr-^'\ (5) 

where .7^ is a constant (Lou 2002; Lou & Fan 2002) and 

Br ^ B^^O . (6) 

For a more general power-law radial variation of the az- 
imuthal magnetic field and those of other related back- 
ground variables, the interested reader is referred to a recent 
analysis of Shen, Liu & Lou (2005). 

In the fluid approximation of a stellar SID, the mass 
conservation, the radial component of the momentum equa- 
tion and the azimuthal component of the momentum equa- 
tion are given below in order, namely 

1 dirT.^u") 1 d{T.'f 
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(8) 



dt dr r2 de de de ' ^ ' 

where u" is the radial component of the bulk flow velocity 
of the stellar SID, j" = rv" is the stellar specific angular 
momentum along the z direction, v" is the azimuthal com- 
ponent of the stellar bulk flow velocity, p is the total grav- 
itational potential, 11" is the vertically integrated pressure, 
and is the surface mass density of the stellar SID. 

In the magnetofiuid approximation for the gaseous 
MSID, the mass conservation, the radial component of the 
momentum equation and the azimuthal component of the 
momentum equation are given below in order, namely 
9E3 19(rE»«») 1 d{T.^f) 
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(10) 



(11) 



(12) 



where is the radial component of the gas bulk flow veloc- 



ity, f 



is the gas specific angular momentum in the z 



direction, is the azimuthal component of the gas bulk flow 
velocity, fl" is the two-dimensional gas pressure, E" is the 
gas surface mass density and Br and Be are the radial and 
azimuthal components of the magnetic field B. The last two 
terms on the right-hand sides of equations Ull^ and H12|l are 
the radial and azimuthal components of the Lorentz force 
due to the coplanar magnetic field. The two sets of fiuid and 
magnetofiuid equations — Q and llOII — 1121 1 are dynam- 
ically coupled by the total gravitational potential p through 
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the Poisson integral 
V{r,e,t) = j dip j 



[(2 + r2 - 2Crcos(^/) - 6)^/^ ' ^^^^ 
The gravitational effect of a massive dark matter halo is not 
included for the moment. The divergence-free condition for 
magnetic field B = {Br, Bg, 0) is simply 



djrBr) dBg 
dr ^ de 







(14) 



and the radial and azimuthal components of the magnetic 
induction equation for gas motions are 

dBr _ 1 a 

dt r I 



-(u^Be-v'Br) , 



(15) 



dBg 9 g 
Ot or 



V^Br). (16) 

Using Poisson integral 1)13^ , one readily derives the following 
expressions for the background surface mass densities 
al{l+Dl) 



2nGr{l + S) 



and 



[al{l + Dl)-Cl/2]5 



(17) 



(18) 



2nGr{l + S) 

where 5 = Eq/Eq is the surface mass density ratio of the 
two dynamically coupled background (M)SIDs and Ca is 
the constant Alfven wave speed in the MSID defined by 



dzBl/iA-K^l) 



(19) 



Apparently, equations 11711 and 1181 requires 

al{l + Dl) ^ aid + Dl) ~ Cl/2 . (20) 

We now introduce two more useful dimensionless parameters 
here. The first parameter /3 = /a^ stands for the square of 
the ratio of the stellar velocity dispersion to the isothermal 
sound speed in the MSID and the second parameter = 
C^/flg stands for the square of the ratio of the Alfven wave 
speed to the isothermal sound speed in the MSID. In late- 
type disc galaxies, the stellar velocity dispersion as is usually 
higher than the gas sound speed ag, we thus focus on the 
case of /3 > 1 (Jog & Solomon 1984a, b; Berlin & Romeo 
1988; Jog 1996; Elmegreen 1995; Lou & Fan 1998b; Lou & 
Shen 2003; Shen & Lou 2003, 2004a, b; Lou & Wu 2005). 

Before going further, we note that a typical disc galaxy 
system involves a massive dark matter halo, a stellar disc 
and a gaseous disc of interstellar medium (ISM) on large 
scales, where the ISM disc is magnetized with the magnetic 
energy density (~ leV/cm'^) being comparable to the en- 
ergy densities of thermal ISM and of relativistic cosmic- 
ray gas (e.g. Lou & Fan 2003). To comprehensively un- 
derstand multi-wavelength observations of large-scale spi- 
ral structures of disc galaxies and to develop potentially 
powerful observational diagnostics (Lou & Fan 2QQQa, b), it 
would be more realistic and necessary to take into account 
of magnetic field effects in a composite magnetized disc-halo 
model.* While there are exceptions, galactic magnetic fields 



The cosmic-ray gas is set aside here merely for the sake of sim- 



typically tend to be coplanar with the disc plane of a spiral 
galaxy on large scales. On smaller scales, regions of closed 
and open magnetic fields are most likely intermingled by 
the solar analogy (e.g. Lou & Wu 2005). As a first step, 
Lou (2002) carried out a coplanar MHD perturbation anal- 
ysis for stationary aligned and unaligned logarithmic spiral 
structures in a single MSID embedded with an azimuthal 
magnetic field and demonstrated that the minimum of the 
ring fragmentation curve in this MSID model is clearly re- 
lated to the generalized MHD Qm parameter (Lou & Fan 
1998a). Since the background MHD rotational equilibrium 
adopted by Lou & Fan (1998a) is not an MSID model, it 
would be more satisfying to justify the statement and in- 
terpretation of Lou (2002) in a dynamically self-consistent 
manner. We shall define a. Qm parameter similar to that of 
Lou (2002) and show that this Qm is equivalent to that of 
Lou & Fan (1998a). 

For a single MSID, we readily derive linear equations 
(by setting relevant parameters for the stellar SID to vanish) 
for coplanar axisymmetric MHD perturbations with Fourier 
harmonic dependence exp(ifcr -|- iojt), where k is the radial 
wavenumber and uj is the angular frequency. In the usual 
WKBJ or tight- winding approximation of fcr 3> 1, we obtain 
the local WKBJ dispersion relation for MHD density waves 
propagating in an MSID in the form of 

= nl + k''{al + C\)-2-KG\k\T.l , (21) 

which is the MHD generalization of the WKBJ dispersion 
relation for coplanar perturbations in an unmagnetized SID. 

To derive an effective Qm parameter for the axisym- 
metric stability (i.e., uj'^ > 0) against MHD perturbations 
with an arbitrary k, the determinant of the right-hand side 
of equation H21|l should be negative for all k. This requires 



Kg{Cl + a 



2x1/2 



> 1 



(22) 



where Qm is an generalized Toomre's Q parameter for a 
single MSID. Inequality 12211 is of the same form as the ef- 
fective Q parameter derived by Lou & Fan (1998a) for fast 
MHD density waves in a rigidly rotating gas disc with a 
different formalism but with the same assumption fcr 3> 1. 
This result establishes the connection found by Lou (2002) 
that the minimum of the ring fragmentation curve is tightly 
associated with the Qm parameter of Lou & Fan (1998a). 

In our low-frequency time-dependent analysis for a 
composite MSID system, we write coplanar axisymmetric 
MHD perturbations with the same harmonic dependence of 
exp(iA;r -f iut). In the WKBJ limit of fcr ^ 1, we obtain 
the local WKBJ dispersion relation for a composite MSID 
system in the form of 



(tJ 



-fc af + 27rG|fc|Eo) 



X [uj"^ - Kg - k^{al + C\) + 27rG|fclEg] 



(23) 



= (27rG|fc|Eo)(27rG|fciEg) . 
As expected, this dispersion relation explicitly shows the 

plicity. Relativistic cosmic-ray electrons gyrating around galactic 
magnetic field give off the observed synchrotron radiation). 
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mutual gravitational coupling between the two SIDs by the 
term on the right-hand side. The first factor on the left-hand 
side is for perturbations in the stellar SID approximated as 
a fluid, while the second factor is for coplanar MHD pertur- 
bations in the gaseous MSID [see expression l|21^ ]. It should 
be emphasized that while dispersion relation H23^ is a lo- 
cal one, the background physical variables are in rotational 
MHD equilibrium in a consistent manner. In particular, Eq 
given by expression 1181 contains the information of back- 
ground magnetic field via C\. 

As noted by Lou & Shen (2003) , the dispersion relation 
derived here for coplanar MHD perturbations in a composite 
(M)SID approach is qualitatively similar to those previously 
obtained by Jog & Solomon (1984), Elmegreen (1995) and 
Jog (1996) in spirit, but one major distinction is that in the 
formulation of our composite (M)SID, the rotation speeds of 
the two SIDs Vs and Vg are different in general. We thus have 
inequality Ks Kg, while in those earlier analyses, Ks = fig 
was presumed a priori following the assumption of Vs = Vg. 

This is conceptually related to the phenomenon of 
asymmetric drift (e.g. Binney & Tremaine 1987) . Physically, 
stellar velocity dispersions mimic a pressure-like force for the 
stellar component, while the thermal ISM gas and magnetic 
pressure forces together act on the magnetized gas compo- 
nent. In the same gravitational potential well determined 
by the total mass distribution, the difference in the stellar 
pressure-like force and the sum of the gaseous and magnetic 
pressure forces would lead to different Va and Vg and thus 
the asymmetric drift. The rare situation of Va and Vg being 
equal may happen under very special circumstances. 



with the two roots and uj'L given by 

1/2 (29) 
± [{Hi + Hif - 4{HiH2 - GiG2)]''^V2 • 

Similar to the proof of Shen & Lou (2003, 2004a), the ojI 
root remains always positive. In contrast, uj^ may become 
negative, leading to axisymmetric MSID instabilities. Sub- 
stitutions of expressions Hi, H2, Gi, G2 and definitions 
1171 1 and 1181 and expression 1201 into equation 1291 1 for 
the minus-sign solution would give uj^ in terms of five di- 
mensionless parameters Da, K = \k\r, 5, (3 and A^, namely 

LuUk) = ^{MK^ + AiK + A,, ^ p"^) , (30) 

where 

A2 = l + l/P + X" IP , (31) 

Ai = -{l + y), (32) 

Ao = 2 + Av+^^^j^ , (33) 

p = BiK^ + B'iK^ + BiK^ + BiK + Bo , (34) 

B4 = [l-il + \^)/P]\ (35) 

Bs = 2(1 + y){5 - 1)[1 - (1 + A^)//3]/(l + 5) , (36) 



3 AXISYMMETRIC STABILITY ANALYSIS 
FOR A COMPOSITE MSID SYSTEM 

We describe below coplanar MHD perturbation analysis for 
the axisymmetric stability of a composite MSID based on a 
low-frequency time-dependent WKBJ approach. Generaliz- 
ing the notations of Shen & Lou (2003) yet with the effect 
of magnetic field included, we here define 

Hi=K'i + k^al - 27rG|fc|Eo , (24) 
H2 = i4 + k''{al + Cl)-2-KG\k\T.l , (25) 
Gi = 27rG|fclEg , (26) 

G2 = 27rG|fc]E^ , (27) 

where Eq and Eq are given by equations 1171 and 1181 . In 
addition to the appearance of C\ in equation 1251 . Eq given 
by expression 1181 also contains the magnetic field effect. 
Dispersion relation 1231 can be cast into the form of 

uj^ - (Hi + H2)lo^ + (H1H2 - G1G2) = , (28) 



B2 = + 2j/ - 3 (8/3 - 4 - + 2A* + 2p\^)/l3^\ , (37) 
Bi = 4(1 + y){l -S)\l- 1//3 + aV(2/3)]/(1 + 5) , (38) 

Bo = 4[1 - 1//3 + AV(2/3)f , (39) 

where y = Z)J. The analysis here parallels that of Shen & 
Lou (2003); the novel magnetic field effect to be explored is 
contained in the dimensionless parameter A'^. 

As a result of the one-to-one correspondence between 
_Dg and dictated by expression 12U1 , it is straightforward 
to derive an equivalent form of expression 13U1 in terms of 
_Dg instead of D^. Mathematical solutions of and 
become unphysical for either < or < or both be- 
ing negative. One can readily show from condition 1201 that 
Dl < for ^ > 1 (Lou & Zou 2004). Therefore, it sufiices 
to require > 0. In the subsequent analysis, we mainly 
use Dl - the square of the rotational Mach number - to 
examine the axisymmetric stability property in a composite 
MSID system. 

By setting A^ = for zero magnetic field in expressions 
H3()|l — H39|l . they all reduce to the corresponding expressions 
for a composite system of two coupled unmagnetized SIDs 
analyzed by Shen & Lou (2003). For scale-free discs more 
general than SIDs, the reader is referred to the work of Syer 
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& Tremaine (1996), Shen & Lou (2004a, b) and Shen, Liu & 
Lou (2005). To derive an effective MHD Qm parameter for 
a composite disc system of one SID and one gaseous MSID, 
we must determine the value of Kmin at which ui^ reaches 
the minimum value. 

3.1 The D^-Criterion in the WKBJ Regime 

Before defining an effective Qm parameter, we first show 
unambiguously the Ds— criterion for axisymmetric stability 
and confirm our earlier interpretations for the marginal sta- 
bility curves in a composite MSID system (Lou & Zou 2004). 

According to solution 1301 1. is a function of K and 
-Dj. By setting ujI. = and assigning values of parameters 
5, (3 and in solution l)30|l . we end up with an equation 
for DI and K. Contours of uj'i. in Di and K with various 
combinations of 5, (3 and are given below to compare 
with our marginal stability curves obtained earlier (Lou & 
Zou 2004; see also Shen & Lou 2003, 2004a, b, Shen, Liu & 
Lou 2005 and Lou & Wu 2005). 

As an example of illustration, we set |m| = 0, 5 — 0.2, 
/3 — 1.5 and A^ = 1 and determine numerically contour 
curves of lj'L in terms of and K as displayed in Fig. 
Physically, the two regions labelled uj^ < in the lower-left 
and upper-right corners are unstable, while the region la- 
belled by uj^ > is stable against axisymmetric coplanar 
MHD perturbations. For comparison, we show the global 
marginal stability curve in a composite MSID system with 
the same parameter values in Fig. |5| (figure 11 of Lou & Zou 
2004), where a is a dimensionless effective radial wavenum- 
ber (see Shu et al. 2000; Lou 2002; Lou & Shen 2003; Lou & 
Zou 2004). In the WKBJ approximation of large K and a, 
the two upper-right solid curves in Figs. Q and |21 show good 
mutual correspondence for the ring fragmentation instabil- 
ity. Thus our previous interpretation for the global station- 
ary axisymmetric MHD perturbation configuration as the 
marginal stability curve is confirmed by the WKBJ analy- 
sis here. In comparison, the correspondence in the small K 
regime is qualitative with apparent deviations; the WKBJ 
approximation works better for a local analysis, while the 
results of Fig. |5| are global and exact without the WKBJ 
approximation. It is clear that this comparison reveals the 
physical nature of the demarcation curves as the axisym- 
metric stability boundaries. 

By an cut contour plot for versus K, the stable range 
of in the WKBJ approximation can be readily identified. 
For example, in Fig. 0with |m| — 0, 5 — 0.2, f3 — 1.5 and 
A^ = 1, the composite MSID system has a stable range from 
= 0.1205 at K ^ 0.5075 to = 6.3428 at A" = 3.5149. 
In the WKBJ approximation, we explore numerically and 
show some qualitative trends of variations for the marginal 
stability curves with parameters S, /3 and A^ in Figs.|21-|3 In 
general, the increase of S and /3 tends to make a composite 
MSID system more vulnerable to instability (compare Figs. 
[HlllandEll, while the increase of the magnetic field strength 
A'^ expands the stable range and reduces the danger of in- 
stabilities (compare Figs. |3 and 2J- 
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Figure 1. A contour plot of cj^ as a function of K and DJ 
with \m\ = 0, 5 = 0.2, /3 = 1.5 and = 1. The two separated 
regions labelled u)'^ < in the lower-left and upper-right corners 
are unstable. The two solid curves mark = 0. 
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Figure 2. The global marginal stability curve of versus ef- 
fective dimensionless radial wavenumber a for \m\ = 0, 5 = 0.2, 
(3 = 1.5 and = 1 [see figure 11 of Lou & Zou (2004)]. 
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Figure 3. A contour plot of lj^ as a function of K and 
with |m| = 0, (5 = 0.2, (3 = 1 and A^ = 0.09. The two separated 
domains labelled by o;^ < are unstable. The two solid curves 
mark u)^ = 0. 
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Figure 4. A contour plot of ui^ as a function of K and with 
\m\ = 0, (5 = 0.2, (3 = 1 and = 3.61. The two domains labelled 
are unstable. The two solid curves mark oj^ = 0. 



Figure 7. The global marginal stability curve of versus ef- 
fective dimensionless radial wavenumber a for |m| = 0, 5 = 0.2, 
13=10 and = 1 [see figure 13 of Lou & Zou (2004)]. 
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Figure 5. A contour plot of uP_ as a function of K and with 
\m\ = 0, (5 = 0.2, (3 = 10 and A^ = 1. The two domains labelled 
are unstable. The two solid curves mark lj^ = 0. 
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Figure 6. A contour plot of iAp_ as a function of K and with 
\m\ = 0, 5 = 1, /3 = 10 and A-^ = 1. The two regions labelled by 
are unstable. The two solid curves mark uP_ = 0. 



The role of coplanar ring magnetic field in stabilizing 
a composite MSID system can be physically understood as 
follows. According to dispersion relation H23|l . we see that, 
in the MSID factor, the gas and magnetic pressure terms are 
explicitly associated with the square of the radial wavenum- 
ber |fc|, while the background surface mass density term is 
linear in |fc|. For the ring fragmentation instability that oc- 
curs at relatively large radial wavenumbers, the two pressure 
terms dominate and the increase of magnetic pressure tends 
to enhance the axisymmetric stability of a composite MSID 
system. For the Jeans collapse instability that occurs at rela- 
tively small wavenumbers, the background gas surface mass 
density term in the MSID factor becomes dominant over the 
two pressure terms. By the rotational MHD radial force bal- 
ance condition IIBII . the background gas surface mass den- 
sity tends to be reduced by the increase of magnetic field 
strength and thus the Jeans collapse instability tends to be 
suppressed. In addition, the right-hand side of dispersion re- 
lation 12311 represents the mutual gravitational coupling in 
the presence of coplanar MHD perturbations. A reduction of 
background gas surface mass density will weaken this cou- 
pling and thus increase the axisymmetric stability. 

For a further comparison, we reproduce the global 
marginal stability results of figure 13 in Lou & Zou (2004) 
here as Fig.|7| which has the same set of parameters as Fig.|K| 
for the local WKBJ solution results. We note again that for 
the unstable region (upper right) of large radial wavenum- 
ber, labelled as the ring fragmentation in Lou & Zou (2004) , 
Figs. 1^ and [7| show very good correspondence as expected. 
As reference. Table 1 contains several lists for the overall 
stable range of Dl with m = Q and different sets of pa- 
rameters 5, fJ and A^, including the results both from here 
using the WKBJ approximation and from those of Lou & 
Zou (2004) for exact global MHD perturbation calculations 
in a composite MSID model. 

As shown above, axisymmetric stability properties of 
a composite MSID system can be qualitatively understood 
using the WKBJ analysis. Nevertheless, the WKBJ approxi- 
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Table 1 . The stable range of for axisymmetric coplanar MHD 
perturbations of any wavelengths for different values of S, (3 and 
A'^ in the WKBJ approximation. The values in parentheses are 
those determined by Lou & Zou (2004) for global marginal sta- 
bility curves and are more accurate, especially for describing the 
Jeans collapse regime involving large radial spatial scales. 



5 


/3 


A2 


lower limit of 


upper limit of 


0.2 


1.5 


1 


0.1205 


(0.7695) 


6.3428 


(6.1554) 


0.2 


1 


0.09 


0.1594 


(0.9063) 


5.9573 


(5.7561) 


0.2 


1 


3.61 


0.1111 


(0.6783) 


7.9494 


(7.7905) 


0.2 


10 


1 


0.1184 


(0.7534) 


4.5499 


(4.4310) 


1 


10 


1 


0.0396 


(0.4063) 


2.2787 


(2.2434) 



mation gradually becomes inaccurate in dealing with smaller 
radial wavenumbers of the Jeans collapse regime. Once the 
physical interpretation has been firmly established, the ex- 
act global marginal perturbation solution procedure of Lou 
& Zou (2004) should be adopted to identify the relevant 
stable range of Dl in a composite (M)SID system. 

3.2 MHD Generalizations of Effective Q 

Parameters in a Composite MSID System 

For the axisymmetric stability of a single disc, the famil- 
iar Q parameter (Safronov 1960; Toomre 1964) provides a 
simple local criterion Q > 1. For a single magnetized disc 
with a coplanar magnetic field, the MHD generalization of 
this local criterion for the axisymmetric stability becomes 
Qm > 1 (Fan & Lou 1996; Lou & Fan 1998a). For a com- 
posite system of two gravitationally coupled discs without 
magnetic field, it would be natural to seek an extension of 
such a local criterion for the axisymmetric stability. Such 
a local criterion does exist but to find a relatively simple 
analytical expression is not as easy (Elmegreen 1995; Jog 
1996; Shen & Lou 2003). It should be emphasized that in 
the analyses of both Elmegreen (1995) and Jog (1996), the 
local background variables are prescribed a priori so that the 
two disc rotation speeds and thus the two epicyclic frequen- 
cies are the same. In contrast, in the analyses of Lou & Shen 
(2003) and Shen & Lou (2003), the equilibrium background 
variables are determined in a dynamically consistent man- 
ner so that the two SID rotation speeds and thus the two 
epicyclic frequencies are allowed to be different in general. 
In reference to the local criteria of Elmegreen (1995) and 
Jog (1996), the global D^-criterion of Lou & Shen (2003) 
for the axisymmetric stability in a composite SID system is 
more straightforward and precise (Shen & Lou 2003). The 
main thrust of this section is to discuss and establish the 
global _Ds— criterion for the axisymmetric stability in a com- 
posite MSID system with a coplanar magnetic field through 
comparisons (Lou & Zou 2004). 

3.2.1 MHD Extension ofQ^i Parameter of Elmegreen 

Parallelling the procedures of Elmegreen (1995) and of Shen 
& Lou (2003), we may define an effective Qe parameter 



as the MHD extension of the Qe' parameter of Elmegreen 
for a composite system of MSIDs. From equation I30L the 
minimum of ui^ is given by 

„2 



1/2, 



where 



K 



^)(Q|-l) , 



(pV2 - A^Kl^^ ~ A,K^,„) 



(40) 



(41) 



and p defined by equation 1)34^ takes on the value at Kmin. 
Although the forms of these mathematical expressions are 
strikingly similar to those of Elmegreen (1995), all relevant 
coefficients and variables contain the effect of magnetic field 
througn A^. By definition 14H of Q% above, it is clear that 
for Q% > 1, the minimum of cot > and thus the com- 
posite MSID system would be stable against axisymmet- 
ric coplanar MHD perturbations of arbitrary radial wave- 
lengths. This generalized MHD parameter Q% corresponds 
to the stable range of where Ds is the rotational Mach 
number of the stellar SID modelled as a fiuid. 

The formidable appearance of cj?. given by equation 
1401 1 would prevent us from deriving a straightforward ana- 
lytical expression of Kmin ■ Nonetheless, instead of minimiz- 
ing oji directly, it is much simpler to determine the critical 
value Kc for K corresponding to the minimum of variable 
W = uj^Lu^. By equation l|28^ . we have 

W = ujIlu1 = H1H2- G1G2 ■ (42) 

For possible extrema of W, the relevant cubic equation that 
Kc should satisfy is 

rfW _ d{HrH2 - G1G2) _ „ 

dK dK ' ^ ' 

or equivalently 

dK^ + aK^ + bK + c = , (44) 
with the four coefficients explicitly defined by 
d = 4 (l + A^) (5+1), 



(45) 



a = -3 {(3S + \^ + l) (y+1), 



(46) 



6 = 2 (5-t-l) {2(3y-2 + X^ + 2(3 + 2y\^ + 2y) , (47) 

c=-{y+l){2/3y + 2PyS-2 + \^ + 2p) . (48) 

For most parameter regimes under consideration, there is 
only one real solution for the cubic equation 14411 . This real 
solution Kc takes the lengthy but straightforward form of 

Kc^{x~ q/2f'^ + {-X - q/2f'^ - a/(3d) , (49) 

where x = {q^ /A + p-^ /27Y''\ p = (b/d) - {a/df/3 and q = 
2(a/d)^ /27- ah/ (3d'^)+c/d. We then use this Kc to estimate 
Kmin and to determine the effective MHD Qe parameter 
as the generalization of Qe'- Relevant curves of MHD Q% 
versus corresponding to different values of 5, /? and 
are displayed in Figures IHl- llOl 
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Figure 8. Several curves of Q|, versus for different /3 values 
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Figure 10. Several curves of MHD Q|, versus for different 
values. Other relevant parameters |m| = 0, 5 = 0.2 and /3 = 1.5 
are fixed. For each curve, the two intersection points at Q? 



bracket the stable range of D j. 
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Figure 9. Several curves of MHD Q|, versus for different 5 
values as indicated. Other relevant parameters |m| = 0, /3 = 10 
and A2 = 1 are fixed. For each MHD Q „ curve, the two intersec- 
tion points at Q|, = 1 bracket the stable range of Dl. This stable 
D| range shrinks as (5 increases. 

By varying parameters (5, (i and A^, we observe several 
trends of variation in the profile of MHD Q\ versus D\ as 
displayed in Figs. |gl-[Tni When /3 increases with fixed (5 and 
A'^ values, or & increases with fixed (i and A^ values, the 
stable range of shrinks in general, while the increase of 
A^ tends to expand the stable range of D\ . Similar variation 
trends have been noticed earlier for the D^— criterion and the 
exact global perturbation solutions in our composite MSID 
model (Lou & Zou 2004). 

By definition 14111 for the MHD Qe parameter, one de- 
termines the stable range of Dl in the WKBJ approxima- 
tion as shown in Figures l51-[Tni For example, given 5 = 0.2, 
/3 = 1.5 and A^ = 1 in Fig. |H1 we have Q| > 1 for 
0.12 < D1 < 6.34; the composite MSID system is thus 
stable against axisymmetric coplanar MHD perturbations 
within this range of D^. More numerical results for the sta- 
ble ranges of D\ for different parameters are summarized in 



Table 2. Approximate stable ranges of determined by the 
criterion > 1 in a composite system of (M)SIDs. The val- 
ues outside the parentheses are given by the MHD parameter 
generalizing that of Elmegreen (1995), while the values in- 
side the parentheses are derived from the MHD Qj parameter 
generalizing that of Jog (1996). These two sets of values cor- 
responding to MHD Qe and MHD Qj parameters are nearly the 
same and the stable ranges are very close to the values given 
by our exact global Ds— criterion using the MHD perturbation 
procedure of Lou & Zou (2004; see also Table IHhereV 



<5 


/3 


A2 


lower limit of 


upper limit of 


0.2 


1.5 


1 


0.12 (0.12) 


6.34 (6.34) 


0.2 


1 


0.09 


0.16 (0.16) 


5.96 (5.96) 


0.2 


1 


3.61 


0.11 (0.11) 


7.95 (7.95) 


0.2 


10 


1 


0.12 (0.12) 


4.55 (4.55) 


1 


10 


1 


0.04 (0.04) 


2.28 (2.28) 



Table|5| By inspection and comparison, the stable D1 ranges 
are almost the same as those derived from our _Ds— criterion 
as expected (see Table 0. In comparison with the stable 
D\ ranges derived from the exact global axisymmetric per- 
turbation solutions of Lou & Zou (2004), we see again that 
the results obtained by the two procedures have very good 
correspondence in the larger wavenumber regime, while in 
the Jeans collapse regime of smaller radial wavenumbers, 
the WKBJ approximation leads to apparent deviations. The 
generalized instability criterion characterized by the MHD Q 
parameter is physically relevant to the MHD ring fragmenta- 
tion instability in a composite disc system of gravitationally 
coupled SID and MSID. 

Alternatively, we may follow a similar procedure of Lou 
& Zou (2004) to show the tight correspondence between the 
marginal stability curve and the effective MHD Q parame- 
ter. For specific values of parameters 5, (i and A'^, we cal- 
culate the minimum of the ring fragmentation curve for 
using the procedure of Lou & Zou (2004). Inserting the re- 
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Table 3. Numerical values for the minima of ring fragmen- 
tation curve (Lou & Zou 2004) and the corresponding values of 
effective MHD Q parameters, including both the MHD pa- 
rameter generalizing Q^/ of Elmegreen (1995) and the MHD Qj 
parameter generalizing Qji of Jog (1996). 
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/9 


A2 




MHD Q| 


MHD Q2 


0.2 


1.5 


1 


6.1554 


1.0216 


1.0223 


0.2 


1 


0.09 


5.7561 


1.0248 


1.0247 


0.2 


1 


3.61 


7.7905 


1.0144 


1.0159 


0.2 


10 


1 


4.4310 


1.0184 


1.0188 


1 


10 


1 


2.2434 


1.0115 


1.0102 



suiting Ds into expression H41^ . we readily obtain the value 
of MHD QIj parameter corresponding to the minimum of 
the ring fragmentation curve. For example, given 5 = 0.2, 
13 = 1.5 and = 1 in figure 11 of Lou & Zou (2004) (or 
Fig. 121 here), the minimum of in the ring fragmentation 
curve is about 6.1554. Using definition I41II . we obtain the 
corresponding MHD Qjj = 1.0216. More numerical results 
are summarized in Table 

In all these cases, the relevant values of MHD are 
fairly close to unity. Therefore the effective MHD Q param- 
eter Qd is indeed closely relevant to the MHD ring frag- 
mentation instability for axisymmetric coplanar MHD per- 
turbations in a composite MSID system. 



3.2.2 MHD Extension ofQji Parameter of Jog 

We have seen that the MHD Qe parameter generalizing the 
Qe' of Elmegreen (1995) is pertinent to the MHD ring frag- 
mentation perturbation in a composite MSID system. How- 
ever, the determination of the MHD Qe parameter may 
occasionally become cumbersome when equation II44II has 
three real roots of K. It would be tedious to identify the 
absolute minimum of cot . Alternative to this procedure. Jog 
(1996) adopted a seminumerical approach to define another 
effective Q parameter for an unmagnetized composite disc 
system, referred to as Qj' here. For comparison, Shen & Lou 
(2003) have used a generalized procedure to derive the Qj 
parameter for a composite SID system with a self-consistent 
background equilibrium allowing two different SID rotation 
speeds but in the absence of magnetic field. We shall intro- 
duce below the MHD Qj parameter for a composite (M)SID 
system with a coplanar azimuthal magnetic field and com- 
pare these results with those of the MHD Qe parameter. 

According to the minus-sign solution 12911 . uj^ becomes 
positive and negative depending on whether H\H2~GiG2 > 
and < 0, respectively. The critical condition for the ax- 
isymmetric stability is thus characterized by HiH2~GiG2 = 
which can be rearranged into the form of 

27rGfcEg 



27rGfcES 



[K2+fe2(a2 + Ci)] + (^i + Pai) 

Note that the form of the expression here is strikingly similar 
to that of Jog (1996) with the magnetic terms being both 
explicit in the expression and implicit in Eg through Cj. 
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Figure 11. Several curves of the MHD Qj versus for different 
13 values. Other relevant parameters |m| = 0, S = 0.2 and = 1. 
For each MHD Qj curve, the two intersection points at Qj = 1 
bracket the stable range of D2. This stable range shrinks as 
(3 increases. 

With magnetic field, we define a new variable by 

27rG'fcEg 



27rGfcEg 



+ 



[Ac2+fc2(a2 + Ci)] (^i + Pai) 

Kf3{l+y)5/{l + 5) 
2[(3{1 +y)-l+ A2/2] + K^1 + A^) 
K(l+y)/{l+5) 



(51) 



i2y + K^) 

where y = and K = \k\r. We then numerically search 
for Kmin at which tot reaches the minimum value. We now 
define an effective MHD Qj parameter such that 



(1 



= J^ = 



j^m»n(l+y)/(l+g) 
Krr^^nf3{l+y)S/il+S) 



(52) 



^ 2[/3(l + y)-l + AV2] + Kl.Jl + A^) 
[see equations (5) and (6) of Jog (1996) and equation (30) 
of Shen & Lou (2003)]. It follows immediately that Qj > 1 
and Q J < 1 correspond to axisymmetric MHD stability and 
instability, respectively. Given Qj > 1, it indeed follows that 
H1H2 — G1G2 > for the luL corresponding to Kmin. That 
is, the minimum of uj'L is positive for arbitrary K and con- 
sequently, the composite (M)SID system is stable against 
axisymmetric coplanar MHD disturbances. The procedure 
of obtaining the MHD Qj parameter can be summarized as 
follows. For a given set oi 5, (3, \^ and D^, one first deter- 
mines the value of Kmin numerically using equation (I30II . By 
inserting this Kmin into equation 1521 , one then obtains the 
numerical value of MHD Qj in a composite MSID system 
for the given set of (5, (3, A^ and D^. 

We have explored relevant parameter regimes of interest 
and revealed several qualitative variation trends of MHD 
Q J parameter as shown in Figures [TT1-[T^ Similar to the 
MHD QI parameter, the range of Dl for MHD Qj > 1 
corresponds to the axisymmetric stability of a composite 
MSID system. Corresponding to MHD > 1, we have 
calculated stable ranges of given several sets of 5, /3 and 
A^ and the detailed results are summarized in Table|5|along 
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Figure 12. Several curves of MHD Qj versus for different <5 
values. Other relevant parameters |m| = 0, /3 = 10 and = 1 
are fixed. For each MHD curve, the two intersection points 
at Qj = 1 bracket the stable range of D^. This stable range 
shrinks as S increases. 
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Figure 13. Several curves of MHD Qj versus for different 
values. Other relevant parameters |?Tt| = 0, (5 = 0.2 and /3 = 1.5 
are fixed. For each Q j curve, the two intersection points at Q j = 
1 bracket the stable range of . This stable range expands 
as A^ increases. 



with those corresponding to the MHD Q% parameter. It is 
apparent that the stable range of MHD given by MHD 
Qj > 1 is almost the same as that given by MHD Q% > 1. 
Meanwhile, in Figures [TT1 - I13I we reveal the same variation 
trends as those obtained by MHD Q% and by exact global 
MHD perturbation procedure of Lou & Zou (2004). When 
(3 increases with fixed S and values (see Fig. Illjl . or S 
increases with fixed f3 and A^ values (see Fig. 112^ , the stable 
range of Dg tends to shrink, while the increase of A^ tends 
to expand the stable range of (see Fig.O. 

Once again, we follow the similar procedure of Lou & 
Zou (2004) to reveal the close relationship between the min- 
imum of the marginal stability curve and the effective 
MHD Q parameter. For specified values of parameters 5, 
P and A^, we compute the minimum of the ring fragmen- 
tation Dj. curve using the exact global MHD perturbation 
procedure of Lou & Zou (2004), insert the resulting Dl into 



expression 15211 and obtain the value of MHD Qj parame- 
ter corresponding to the minimum of the ring fragmentation 
curve. For example, given 5 = 0.2, (3 = 1.5 and A^ = 1 in 
figure 11 of Lou & Zou (2004) (or Figure |21 here), the min- 
imum of in the ring fragmentation curve is ~ 6.1554. 
Using definition l|52|l . we obtain the corresponding MHD 
Qj — 1.02223. More detailed numerical results are summa- 
rized in Table 13 for reference. 

By these numerical experiments, we demonstrate that 
the values of MHD Qj and MHD Q% corresponding to the 
minima of Dl ring fragmentation curves are nearly the same, 
with the relevant values of MHD Qj being very close to 
unity. Therefore, the effective MHD Qj parameter is also 
pertinent to the MHD ring fragmentation instability for 
axisymmetric coplanar MHD perturbations in a composite 
MSID system. 

In comparison to the determination of MHD Q%, the 
search for the MHD Qj value requires numerical explo- 
rations for each given Dl. The major advantage is that the 
definition l|52^ for MHD Qj remains valid for the entire pa- 
rameter regime and avoids improper situations of unusual 
sets oi 5, 13 and A^. 



3.3 A Composite Partial MSID System 

In most disc or spiral galaxies, there are overwhelming ob- 
servational evidence for the existence of massive dark matter 
haloes in general. To include the large-scale gravitational ef- 
fect of a massive dark matter halo, we add a gravitational 
potential <E> term associated with the dark matter halo in 
our basic MHD equations Q , @ , lITTl and llT^ , where <I> is 
presumed to be axisymmetric for simplicity and for the lack 
of information. Based on A'^— body numerical simulations for 
galaxy formation, typical velocity dispersions of dark matter 
'particles' are fairly high (more than a few hundred kilome- 
ters per second). Hence, another major simplification of our 
analysis is to ignore perturbation responses of the massive 
dark matter halo to coplanar MHD perturbations in a com- 
posite MSID system (e.g., Syer & Tremaine 1996; Shu et al. 
2000; Lou 2002; Lou & Fan 2002; Lou & Shen 2003; Shen 
& Lou 2004a, b). As before, we introduce a dimensionless 
ratio F = (fi/{(f + for the fraction of the disc potential 
relative to the total potential all in a background equilib- 
rium state (e.g., Syer & Tremaine 1996; Shu et al. 2000; 
Lou 2002; Lou & Shen 2003; Lou & Zou 2004). The back- 
ground rotational MHD equilibrium of a composite MSID 
is thus strongly modified by this additional <1> term. As be- 
fore, we write Qs ~ asDs/r, Qg — agDg/r, Ks = \/2Q,s and 
Kg — \/2Q,g. It follows from the radial force balance in the 
(M)SID system that the background surface mass densities 
now become 
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Figure 14. An uP_ contour plot as a function of K and with 
m = 0, <5 = 0.2, p = 1.5, = 1 and F = 0.1. The region labelled 
by u)'^ < is unstable and the Jeans collapse regime disappears 
completely. In comparison with Fig. it is apparent that the 
stable range of is enlarged as F decreases. 

where < F < 1 for a partial composite MSID and F = 1 
for a full composite MSID that has been studied in detail 
in the previous subsections. Performing the standard MHD 
perturbation analysis, we linearize dependent physical vari- 
ables but ignore dynamical feedbacks from the massive dark 
matter halo to coplanar MHD perturbations in the MSID 
system. In the WKBJ approximation, it is then straightfor- 
ward to derive a strikingly similar dispersion relation in the 
form of equation ^2?,\ but with modified background equilib- 
rium properties 1)53^ and H54|l . Following the same procedure 
of Ds— criterion analysis described in section 4.1, we obtain 
the u)'L contour plot as a function of stellar rotational Mach 
number Di and radial wavenumber K = \k\r with the po- 
tential ratio F as an additional parameter. Typical results 
are displayed in Fig. 1141 as an example of illustration. 

The example of Fig.Qwith F = 1 corresponds to a full 
composite MSID system where no background dark matter 
halo is involved, as has been studied in subsection 3.1. As 
F becomes less than 1 corresponding to an increase of the 
potential fraction of the massive dark matter halo, the stable 
range of becomes enlarged as clearly shown in Figure IT^ 
for the case oi F — 0.1. From these oji < contour plots, it 
is apparent that the introduction of a massive dark matter 
halo tends to stabilize a composite MSID system as expected 
(Ostriker & Peebles 1973; Binney & Tremaine 1987; Lou 
& Shen 2003). For late- type spiral galaxies, one may take 
F = Q.l or smaller. Such a composite partial MSID system 
are stable against axisymmetric coplanar MHD disturbances 
in a wide range of D^. 

3.4 Quantitative Estimates and 
Galactic Applications 

For our numerical examples, the disc mass density ratio pa- 
rameter 5 = Eq/Eq has been taken to be 0.2 and 0.02 (see 
Fig. ll2ll . For late-type spiral galaxies, this ratio ranges from 
0.05 to 0.15. For relatively young and gas-rich spiral galax- 



ies, this ratio S can reach 0.2 and even higher. In nearby 
spiral galaxies, the strength of magnetic field is typically in- 
ferred to be a few to lOfiG; by the equipartition argument, 
magnetic field strength may reach a few tens of in cir- 
cumnuclear regions and within towards the centre (e.g., Lou 
et al. 2001). We shall take the ratio A = CA/ag to be of the 
order of 1. As estimates, we have taken /3 = (as/ag)^ to be 
1, 1.5, 10, and 30 in Figure ITTl For spiral galaxies, the ratio 
as/ttg can be of the order of or greater than 5 or 6. With 
these estimates in Figs. 171- 1121 we see clearly that without a 
massive dark matter halo, a typical disc galaxy system with 
a sufficiently fast rotation (e.g., a Vs ~ 150 — 250 km s~^) 
would suffer the ring-fragmentation instability. Although the 
presence of magnetic field offers the stabilizing effect (see 
Fig. against the ring-fragmentation instability, the devel- 
opment of such instability would be unavoidable for typi- 
cally inferred galactic magnetic field strenghs. By Fig. 1141 
the inclusion of a massive dark matter halo holds the key to 
prevent such ring-fragmentation instability. 

In the observational studies of Kennicutt (1989) on 
global star formation rates in spiral galaxies, the theoret- 
ical rationale is therefore not sufficiently strong. Firstly, the 
application of the Toomore stability criterion for a single 
disc is too simple to be physically sensible. Secondly, for a 
composite disc system without magnetic field, the general- 
ized criteria (Elmegreen 1995; Jog 1996; Lou & Fan 1998) 
for the Toomre instability cannot be readily applied to a real 
galactic disc. Thirdly, even for a composite disc system with 
a coplanar magnetic field which shows clear stabilizing ef- 
fects, the MHD ring-fragmentation would occur for typically 
inferred parameters of a spiral galaxy without a sufficiently 
massive dark matter halo. Finally, our conclusion is there- 
fore that in relating Toomre-type instabilities with global 
star formation rates in spiral galaxies, one should conceive 
new physical rationales by incorporating the dynamical in- 
terplay between a massive dark matter halo and a magne- 
tized composite disc system (Wang & Silk 1994; Silk 1997; 
Lou & Fan 2002a, b). 



4 SUMMARY AND DISCUSSION 

In this paper, we examined the axisymmetric MHD linear 
stability properties of a composite system consisting of a 
stellar SID and a gaseous MSID coupled by the mutual 
gravity and a massive dark matter halo, using the WKBJ 
approximation. Our main purpose is to confirm the physical 
interpretation for the global marginal stability curve (Lou 
& Zou 2004) and to establish the MHD generalization of 
the Q parameter (Safronov 1960; Toomre 1964) for a com- 
posite MSID system in reference to the earlier work (e.g., 
Elmegreen 1995; Jog 1996; Lou & Fan 1998a, b; Lou 2002; 
Lou & Shen 2003; Lou & Zou 2004). We now summarize the 
main theoretical results below. 

We have recently constructed exact global solutions for 
stationary coplanar MHD perturbations in a composite sys- 
tem of a stellar SID and a gaseous MSID for both aligned 
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and unaligned logarithmic spiral cases (Lou & Zou 2004). 
Lou & Zou (2004) have extended the analyses of Shu et al. 
(2000) on an isopcdically magnetized SID, of Lou (2002) 
on a single coplanarly magnetized SID, and of Lou & Shen 
(2003) and Shen & Lou (2003) on an unmagnetized compos- 
ite SID system. In a broader perspective, a composite MSID 
system is only a special case belonging to a more general 
class of composite scale-free magnetized disc systems (Syer 
& Tremaine 1996; Shen & Lou 2004a, b; Shen et al. 2005). 
In analogy of Shu et al. (2000), Lou & Shen (2003) and Shen 
& Lou (2003), Lou & Zou (2004) naturally interpreted 
marginal curves for stationary axisymmetric coplanar MHD 
perturbations with radial propagations as the marginal sta- 
bility curves. Based on the low-frequency time-dependent 
WKBJ analysis here (Shen & Lou 2003, 2004a), we estab- 
lish the physical scenario for the presence of the two unstable 
regimes referred to as the 'MHD ring fragmentation insta^ 
bility' and the 'MHD Jeans collapse' in a composite MSID 
system. Consequently, it is intuitively appealing and phys- 
ically reliable to apply our exact global Dg— criterion for a 
composite MSID system in order to examine its axisymmet- 
ric stability and obtain the stable range of Di, where Ds is 
the rotational Mach number of the stellar SID. 

In the WKBJ approximation, we relate our MHD 
Ds— criterion to the two effective Q parameters extended to 
the MHD regime, namely, the MHD Qb— criterion general- 
izing that of Elmegreen (1995) and the MHD Qj— criterion 
generalizing that of Jog (1996). However, our procedures 
differ from those of Elmegreen (1995) and of Jog (1996), be- 
cause our MHD background of rotational equilibrium is dy- 
namically self-consistent with Ks ^ Kg in general. We show 
that MHD generalizations of both Qe— and Qj— criteria 
lead to nearly the same stable range for extended to the 
MHD realm. This confirms the close relation between our 
MHD Da— criterion and the MHD Qe— and Qj— criteria. 
Complementarily, we show that the values of the MHD Qe 
and Qj corresponding to the minima of the DJ ring frag- 
mentation curves in the exact global MHD perturbation pro- 
cedure of Lou & Zou (2004) are all close to unity. Our in- 
terpretation of the ax;isymmetric marginal stability curve 
as the demarcation of stable and unstable regimes is phys- 
ically sensible, and the MHD Qe and Qj parameters arc 
associated with the MHD ring fragmentation instability in 
a composite MSID. 

Finally, we have shown the axisymmetric MHD stability 
property of a composite partial MSID by including the grav- 
itational effect from an axisymmetric massive dark matter 
halo. It is apparent that the dark matter halo has a strong 
stabilizing effect for a composite MSID system. 

In addition to theoretical interest of disc instabili- 
ties for forming large-scale galactic structures (n.b., non- 
axisymmetric ones are not studied here), there has been a 
keen desire to somehow relate such instabilities to global star 
formation rates in disc galaxies and their evolution (e.g.. 
Jog & Solomon 1984a, b; Kennicutt 1989; Wang & Silk 
1994; Silk 1997). The overall chain of star formation pro- 
cesses from large-scale disc instabilities, to giant molecular 



clouds, to cloud collapses, to clusters of stars, to disc ac- 
cretion onto individual stars and so on is quite complicated 
and involves many scales of different orders of magnitudes. 
Conceptually, large-scale axisymmetric ring structures in a 
disc must be further broken down non-axisymmetrically into 
smaller pieces in order to initiate this conceived chain of col- 
lapses. While various stages of this 'chain' have been inten- 
sively studied separately, the ultimate relation or connection 
between the largo-scale axisymmetric instabilities and the 
global star formation rate remains unclear and needs to be 
established (e.g. Elmegreen 1995; Lou & Bian 2005). 

Observationally, Kennicutt (1989) attempted to infer 
an empirical relation between the Q parameter of the stellar 
disc alone and the global star formation rate. Wang & Silk 
(1994) pursued a similar idea with an estimate of Q param- 
eter for a composite disc system of two fluid discs; however, 
their approximation for Q parameter may be off too much 
under various relevant situations (Lou & Fan 1998b, 2000a). 
Should this line of reasoning indeed contain an element of 
truth for addressing the issue of the global star formation 
rate, then the Q parameter adopted should really correspond 
to that of a composite disc system with a magnetized gas 
disc component and in the presence of a massive dark mat- 
ter halo. The basic physical reason behind this suggestion 
is that stars form directly in the magnetized gas disc under 
the joint gravitational influence of the dark matter halo, the 
stellar disc and the magnetized gas disc itself. If this line of 
reasoning does indeed make physical sense, then an interest- 
ing possibility arises. That is, the dark matter halo may play 
an important role of regulating global star formation rates 
in disc galaxies and thus galactic evolution. For example, if 
the mass of a dark matter halo is very much greater than 
the mass of a composite disc system, then star formation 
activities become weaker. On the other hand, if the dark 
matter halo is not sufficiently massive, then the disc system 
rapidly evolves into a bar system. It is also possible that the 
dark matter halo is only marginal to maintain a stability 
of a composite disc. In this case, the global star formation 
activities in the disc system proceed in a regulated manner. 

Our MHD model analysis in this paper, highly idealized 
in many ways, does contain several requisite elements for es- 
tablishing an MHD generalization of the Q parameter and 
the corresponding criterion for axisymmetric stability or in- 
stability. Observationally, it would be extremely interesting 
to examine the relation between the MHD Qm parameter 
in galactic systems and global star formation rates. This is 
not expected to be a trivial exercise given various sources of 
uncertainties. 

By presuming that such axisymmetric disc instabilities 
characterized by either D^ or Qm parameters might some- 
how hint at or connect to the global star formation rate, 
there are a few model problems similar to the current one 
that can be explored further. For example, the models of 
Shen & Lou (2004b) and Shen, Liu & Lou (2005) can be 
combined to construct a composite disc system consisting of 
two scale-free discs with the gaseous one being coplanarly 
magnetized in the presence of a dark matter halo. Likewise, 
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the work of Lou & Wu (2005) can be generalized to two 
coupled scale-free disc with the gaseous one being isopedi- 
cally magnetized in the presence of dark matter halo (Lou 
& Wu 2006 in preparation). The real situation is more com- 
plicated. We hope these analyses may offer certain hints and 
insights for different aspects of the problem (Lou & Bai 2005 
in preparation). 
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